Gapped solitons and periodic excitations in strongly coupled BEC 
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It is found that localized solitons in the strongly coupled cigar shaped Bose-Einstein condensate 
form two distinct classes. The one without a background is an asymptotically vanishing, localized 
soliton, having a wave-number, which has a lower bound in magnitude. Periodic soliton trains exist 
only in the presence of a background, where the localized soliton has a VF-type density profile. 
This soliton is well suited for trapping of neutral atoms and is found to be stable under Vakhitov- 
Kolokolov criterion, as well as numerical evolution. We identify an insulating phase of this system 
in the presence of an optical lattice. It is demonstrated that the VF-type density profile can be 
precisely controlled through trap dynamics. 
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I. INTRODUCTION 



Cigar shaped Bose-Einstein condensate (BEC) pro- 
vides an ideal venue for realizing different types of soli- 
tonic excitations in one dimension. Dark 0, bright 
[H-Q , grey 0] solitons and periodic soliton trains Q have 
been experimentally observed. In the weak coupling 
regime, the Gross-Pitaevskii (GP) equation, describing 
the mean field BEC dynamics, reduces to the well stud- 
ied non-linear Schrodinger equation (NLSE), possessing 
soliton solutions. The fact that NLSE is an integrable 
system explains the observed soliton dynamics 
In one dimension, the possibility of using these localized 
excitations for producing atom laser and other applica- 
tions [r2T - [l4j , requiring macroscopic coherence, has made 
cigar-shaped BEC an area of significant current interest 
[T5l - [l9j . The fact that scatte ring length can be controlled 
through Feshbach resonance [20( , has given access to both 
weak and strong coupling sectors, as well as attractive 
and repulsive domains. The regulation of the transverse 
trap frequency uj± can also be used to control the cou- 
pling in the cigar-shaped BEC, as has been demonstrated 
experimentally in the generation of the Faraday modes 



The strongly coupled cigar shaped BEC, in the repul- 
sive domain, is characterized by a quadratic non-linearity, 
following the Thomas- Fermi approximation [22l . [23j . In 
contrast to the weak coupling domain, the solitons in the 
strong coupling sector are not well studied. The grey soli- 
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ton dynamics has been numerically investigated, which 
yielded a structure similar to that of the weak coupling 
regime [24]. Here, we study dark solitons, as well as soli- 
ton trains in this non-linear system and find significant 
differences, with the weak coupling sector. We restrict 
ourselves to the repulsive domain, keeping in mind the 
fact that strong coupling BEC will be prone to instability 
in the attractive sector. 

It is observed that unlike the case of NLSE, the soli- 
tons exist in two distinct classes. Starting from a gen- 
eral ansatz solution with non- vanishing background, it 
is found that there can be two domains of localized so- 
lutions. The one with background is a W-type soliton, 
which can exist for all the values of the momentum of the 
envelope profile. The other solution is an asymptotically 
vanishing soliton, without a background, which requires 
a finite momentum to get excited. Both the solutions are 
shown to be stable under the Vakhitov-Kolokolov (VK) 
criterion (25|. It is observed that dependence of the ef- 
fective chemical potential on the profile width separates 
these two solution domains and also ensures their sta- 
bility. The W-type soliton is shown to be dynamically 
stable under the numerical Crank Nicholson finite differ- 
ence method. Interestingly, periodic cnoidal waves can 
only exist with a background. Using a general Pade type 
ansatz, we identify more general solutions. It is shown 
that the background must be real and non-vanishing for 
the W-type solitons to exist. Since W-type soliton is well- 
suited for trapping of atoms, we illustrate the procedure 
for coherent control of this structure, through scattering 
length and trap frequency. 

The paper is organized as follows. In Sec. II, we obtain 
exact localized and periodic solutions of the strongly cou- 
pled GP equation in one dimension, and compare them 
with the weak coupling sector. In Sec. Ill, a more general 
ansatz is employed to identify some of the non-linear ex- 
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citations, unavailable through standard procedure. It is 
found that a unique localized solution can exist only with 
a real background. We also explore the structure of the 
solutions in the presence of an optical lattice and identify 
an insulating phase. Subsequently in Sec. IV, their coher- 
ent control is analytically demonstrated, where both the 
amplitude and width of the IF-type soliton can be ex- 
actly regulated. We conclude in the final section with a 
number of directions for future work. 



II. EXACT SOLITONS OF THE STRONGLY 
COUPLED BEC 

The three dimensional GP equation describing the dy- 
namics of BEC in a cylindrical harmonic trap, V = 
hMu>±(x 2 +y 2 ), is given by 



ih 



0#(r, t) 



— V» + tW r,t)|» 



V 



*(r,t), (1) 



where Uq = 47rfi, 2 a/M, is the effective two body inter- 
action, a being the scattering length. In the quasi-one 
dimensional limit, the condensate wave-function can be 
factorized: 



■*(v,t)=f(z,t)G(x,y,a), 



(2) 



where <r{z) is the local particle density. G(x, y, a) is the 
normalized equilibrium wave function for the transverse 
motion: 



a(z)= / dxdy\^(x,y,z)\ 2 = \f(z,t)f 



(3) 



In the repulsive, strong coupling limit (jiqUq >> Hlu±, 
where ui± is the confining trap frequency), one uses the 
Thomas-Fermi approximation for the transverse profile, 
leading to the condensate equation [22|, 



ih—f(z,t) 



^|1 + 2^V/ 2 (|/M)|-^ 



/(*,*)• 



(4) 



Here, do is the equilibrium density of the atoms far 
away from the axis. The non-linear excitations of this 
system can be probed through an ansatz solution, 



/(*,«) 



Q i(kz— u)t) 



p(0, 



(5) 



with a fast moving component and slowly varying en- 
slope profile 
with satisfying, 



velope profile Here, £ = a(z — vt) and v = jjk, 



leading to, 
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(9) 



along with a relation between the width a, and the 
effective chemical potential e, 



a 2 p" 



99 



ep = 0, 



(6) 



where g = -4Mu i a 1 / 2 /K and e = 2Mu/h + 
4Mu} ± {a a)^ 2 /h- k 2 . 

It is straightforward to check that the following ansatz 
solution, 



p{£) = A + Bcn 2 (tm), 



(7) 



solves Eq.6, where cn(£, to) is the cnoidal function, 
with to being the modulus parameter (0 < m < 1) [26j j . 
and A, B are constant parameters to be determined. We 
note that cn(£, 0) = cos(£) and cn(£, 1) = sech(^). A 
serves here as the background. On substitution, we ar- 
rive at the consistency conditions, 



2Ba 2 (m-l)A + gA 2 + eA = 0, 
-4a 2 (l~2m)+2gA + e = 0, 
— 6to« 2 + gB = 0, 



(8) 



e 2 = 16a 4 (to 2 — to + l) 



(10) 



The sign of e is not fixed, as both the roots of the above 
equation are allowed. The general periodic solution can 
then be written as, 



2a 2 



p(£) = (1 — 2m) =F \J to 2 — to + 1 + 3TOcn 2 (£, m) 

9 L 



representing a soliton train. Consideration of the spe- 
cial case to = 1 yields localized solutions. It corresponds 
to two specific values, e = ±4a 2 . The positive root re- 
quires presence of the background A, and the envelope 
takes the form, 



p(0 = — 

9 



1 - -sech 2 (0 



(11) 
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representing a localized W-type soliton. 

The Vakhitov-Kolokolov criterion points out that the 
integral N(e) — J |p(£)| 2 d£, when varied with respect to 
the effective chemical potential e, indicates the stability 
of the solution p(£). In the present case, 

dN(e) _ 6e 
de g 2 ' 

requiring that e > for the stability of the solution, 
which is consistent with e = 4a 2 . 

For this case we find, 

k 2 >2^-\e\, (12) 

setting a lower limit 5jyi for to, if it is positive, in order 
to generate the W-type soliton. Such bound is not there 
if the driving frequency is negative. 

For the negative root of Eq.10, with m = 1, the back- 
ground vanishes, and we get, 

P(0 = ~sech 2 (0, (13) 
which, under the Vakhitov-Kolokolov criterion, yields, 

dN(e) e 



and hence is stable as e < 0. In this case, k 2 > 
|e| + 2-^, indicating that a finite wave number is needed 
to excite the solution for the positive frequency case. 
Such condition does not arise if the driving frequency is 
negative. This type of velocity restricted solitons have 
been identified in higher order non-linear Schrodinger 
equation relevant for optical fiber pulses in the femtosec- 
ond domain [27j . 

We have numerically evolved the W-type solution, us- 
ing the Crank Nicholson finite difference method, which 
is unconditionally stable. In this analysis, the initial 
profile has been taken as ip(z,t = 0) = ii(z t t = 0) + e, 
where e is a function, which assumes a random value 
at each point. The analysis was carried out with 
dz = 0.0003, dt = 0.0001 for 1000 cycles. Fig. Q] shows 
that, the H^-type soliton remains unchanged with its 
form showing minor perturbation, where e is taken to 
be 10 percent of the peak value of ip. The minima 
positions and the width remained unaltered. We have 
also checked that the evolution was unitary conserving 
the number of particles, upto second order in dt. 

III. COMPLEX ENVELOPE BLOCH SOLITONS 

A general ansatz for identifying complex envelope 
Bloch type solitons can be written as, 

p{0 = {A + iB) + (C + iD)cn 2 (£,m), (14) 

where A,B,C,D are real. The corresponding consis- 
tency conditions lead to, 



g(A + iB)y / A 2 + B 2 — e(A + iB) + 2a 2 (C + iD){m - 1) = 0, 

[g 2 {A 2 + B 2 ) - e 2 ] (A + iB)(C + %D) + g 2 {AC + BD)(A + iB) 2 - 8a 4 (C + iD) 2 (m - 1)(1 - 2m) 

-2ea 2 [(C + iD) 2 (m-l) + 2(A + iB)(l-2m)(C + iD)] = 0, 

g 2 [(C 2 + D 2 ){A + iB) 2 + {A 2 + B 2 ){C + iD) 2 + 4{AC + BD)(A + iB){C + iD)] - e 2 (C + iD) 2 

-4ea 2 {C + iD) [3 (A + iB)m2(C + iD)(l - 2m)+] - 8a 4 (C + iD) 2 [2(1 - 2m) 2 + 3m(m - 1)] = 0, 

g 2 [{AC + BD) + {C -iD){A + iB)]-Qea 2 m-2Aa A m{l~2m) = 0, 

g 2 (C 2 + D 2 )-36a 4 m 2 = 0. (15) 

I 



Solutions of these consistency conditions can be found through a tedious, but straightforward calculation; 
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- B 2 , D = -B, 
PJ P 



B 2 = [p{g^ + tf -Ara A {m- l) 2 ] / [(.g^ + e) 2 - 4a 4 (m - l) 2 {r/pf 



(16) 



where, 



P 

q 



\e 2 + 6(1 - 2m W + 4a 4 (llm 2 - 11m + 2) 

3g 



-, r 



3a m 



4a 2 (l - 2m) + e] ; 



(17) 



r 



and, e 2 = 16a . In case of localized solitons, i.e., for 
m = 1, we arrive at the interesting condition that A — 
= C, leading to the previously found W-type soliton. 
However, periodic cnoidal wave solutions can exist with 
complex background parameters. 

We now explore more general solution space, through 
a Pade type ansatz (28l |. 



A, B, and C are real parameters to be determined 
from the consistency conditions, with AC — B ^ 0. It 
can be seen that these conditions are four in number, 
indicating the constrained nature of the general solution. 
The following consistency conditions, with a 2 = /3 can 
be deduced: 



p(0 



A + Bf(Q 



(18) 



A 2 g + 2BC {m- 1) f3 + A (-2 C 2 (m - l)/3 + e) = 0, 

A 2 Cg + B((2m-l)f3 + e)+A(2Bg + C(P~2m(3 + 2e)) = 0, 

B 2 g + AC 2 (—/3 + 2mf3 + e) + BC {2Ag + f3 — 2m/3 + 2e) = 0, 

B 2 Cg- 2Bm(3 + 2ACm(3 + BC 2 e = 0. (19) 



We concentrate on the localized solutions, because of 
its physical interest: 

e f l-2sech(2Q \ 
P ^^~-g{l + sech(20)- (20) 

Like the previous case, the width and the amplitude of 
the solution are coupled. It is easy to see that the profile 
can be rewritten as, 

p^)=-i(l-^sech 2 (0), (21) 

which is identical to the W-type soliton. 

The above solution can be identified as the unique 
separatrix in the phase space of the solutions of Eq.(6), 
separating the periodic solutions with closed orbits from 
the unbounded ones, represented by open orbits (29j . 



The presence of periodic solutions motivates us to ex- 
plore the nature of the solution in the presence of an 
optical lattice, V(x) = Vq cos 2 x, which necessitates the 
solutions to possess sinusoidal character. For the purpose 
of comparision j30| . the 1-D GP equation is normalized 
in the form, 

id t Tp(z,t) = (~d zz + g\ip(z,t)\ + V(x)- fi)iP(z,t), (22) 

where, g is the normalized two-body interaction. The 
solution is found to be of the form, ip(z,t) — (a + 
&cos 2 z)e~ lut . The parameter values for this insulating 
phase is given by, a = Vo/2g and b — —Vo/g, with 
u = 1 + 1/2 Vq — fx. In contrast to the weak coupling 
case [3l|, where, analytic solutions have been obtained 
for both superfluid and insulating phases, for the present 
case, we have been able to identify only an insulating 
phase, analytically. It is worth observing that in this 
case, the competition between the lattice potential and 
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FIG. 1: (Color Online) Direct numerical evolution of type 
soliton with dz = 0.0003, dt = 0.0001 for 1000 iterations. 



tl>(g,t) = B(t)F(C)e[ 4 *^+2 G W], 



(24) 



where T = A(t){z - l(t)}, l(t) = J*v(t')dt' and 
G(t) — L n(t )dt . The phase has been taken in the 
form $(z, t) = a(t) + b(t)z — ^c(t)z 2 7 where a(t) is a z 

independent phase term: a(t) = ao + f A 2 (t')dt' . 
The solutions are necessarily chirped in time and space, 
with time varying amplitude and width. b(t) is a time 
dependent momentum and c(t) balances the oscillator, 
leading to the Riccati equation: 



the nonlinearity yields sinusoidal solutions, whereas, for 
the localized soliton solutions, the nonlinearity and the 
dispersion are responsible for the existence of the solu- 
tions. In the repulsive domain, the above solutions exist 
both for positive and negative values of Vb- The cor- 
responding energy (in dimensionless unit) is found to be 
E = 4<i7rV 2 + (i7rV 3 -2(i7r/iV r 2 , where d = l/8g 2 . We note 
that the contribution from the interaction term does not 
explicitly contribute to the energy, although the coupling 
parameter appears in E through the solutions. The av- 

v 2 

erage atom number density, v = g^-, is a constant and 
can be controlled by tuning the lattice potential and the 
scattering length. 

IV. COHERENT CONTROL OF THE SOLITONS 

We now investigate the effect of time dependent non- 
linearity and gain or loss on the W-type soliton profile. 
The two points of this soliton, where the order parame- 
ter vanishes, may be useful for trapping of neutral atoms. 
The barrier height and the locations of the minima can 
be controlled, by changing the frequency and the scatter- 
ing length, which can be manipulated through Feshbach 
resonance [32l . l33j . Keeping this in mind, in the folloc- 
ing, we study the control of this localized excitation in 
the following section. For the sake of comparison with 
the weak coupling case [33 - l36| . the GP equation is cast 
in the form 

idtf = -\d 2 z J + 7 (i)|V# + \M{t)z\ + (23) 
Here ip, t and z have been scaled, respectively by 

1/2 

a B , uj±_ and l/aj_, making them dimensionless. 7 = 
2(a(t)/as) 1 / 2 is time dependent non-linearity coeffi- 
cient, controllable through Feshbach resonance; M(t) = 
uj 2 (t)/uj 2 _ is related to the axial trap frequency, which 
can be made time dependent. n{t) = rj(t)/huj± is time 
dependent loss/gain. The oscillator length in the trans- 
verse direction is defined as a± = (fi/Mwi) 1 / 2 and is 
the Bohr radius. We consider the following ansatz solu- 
tion for: 



c t -c 2 (t) = M(t). (25) 

The above can be cast as the familiar Schrodinger eigen 
value equation: 

-0"(t)-M(t)0(i)=O, (26) 

via a change of variable: c(t) = — dln ^> _ The constant 
part of M(t) acts as the eigen value. A number of vari- 
ations in the oscillator frequency can be analytically in- 
corporated from solvable quantum mechanical problems. 
The oscillator can also be made expulsive. The location 
of the condensate profile satisfies, dl(t)/dt — c(t)l(t) = 
b(t). Other parameters are obtained, using the following 
consistency conditions: B(i) = B exp(l/2 f*c(t')dt), 
7 (i) = _i 7o A(t)B(t)exp(~G(t)/2), and A(t) = B 2 (t) = 
b(t). The real part of the GP equation yields, 

F"(T)+gF 2 (T) + eF(T) = 0, (27) 

which takes the form of Eq. |6l with g = 

- 2 A(t)B(t) eG{t)/2 > e = - X and T = A (t)[z - l{t)\. One 
can find various singular and non-singular solutions for 
F(T), using the previous procedure, provided g is con- 
stant. We now illustrate two specific cases of interest. 

I. First, we consider the condition M(t) = Mo 2 , with 
a enstant oscillator frequency. This yields a periodic so- 
lution, with period 7r/2M ; 

1>(z,t) - -- y/Aosec(M t) [1 - ^-sech 2 (T/2)] 
9 ^ 

X e i$(z,t)+iG(t) (28) 

The soliton can be compressed and accelerated through 
the time dependence of coupling. The presence of 
sec(M t) in the amplitude and width of the soliton pro- 
file, leads to its compression and amplification, as time 
increases. 

II. We now turn our attention to the solution in an 
expulsive potential (M(t) = — M5 2 ), where the W type 
density profile is given by, 
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FIG. 2: (Color Online) Temporal behavior of W type soliton, 
without gain/loss and with Mq — 0.5, Ao — .01, A = —1 and 
70 = 10. 

ip(z,t) = -- y/A sech(M t) [1 - ^sech 2 (T/2)} 

x e m*,t). (29) 

This profile has a transient character and the behavior 
in time for the expulsive potential is very different from 
the regular case. 

The fact that these minima locations can be controlled, 
including the barrier height, makes these solutions poten- 
tially attractive, for trapping of neutral atoms and their 
manipulation. The temporal evolution of the density pro- 



file is shown in Fig. 2. One observes gradual reduction in 
amplitude as time progresses. 



V. CONCLUSION 

In conclusion, the cigar shaped BEC in the strong cou- 
pling sector leads to two different types of stable localized 
solitons, absent in the weak coupling regime. Localized 
solitons, with no background, require a finite momentum 
to exist for positive driving frequency. No such restric- 
tions are there for the soliton trains which can propagate 
only on a background. The solutions are found to be 
stable under VK criterion. It is shown that the local- 
ized soliton can be effectively controlled through, scat- 
tering length and trap frequencies, which makes it useful 
for trapping of atoms. It is interesting to note that the 
asymptotically vanishing localized solution is also veloc- 
ity restricted. This type of solutions exist in higher order 
non-linear equations, relevant for femto-second pulses in 
optical fiber [27|. We hope that the velocity restricted 
solitons, so far unobserved in the optical fiber system, 
may find experimental verification in cold atoms. It 
will be of interest to find exact solutions for the non- 
polynomial mean field equation relevant for cigar shaped 
BEC [23[ and compare the transition between strong and 
weak coupling sectors. It will be of deep interest also to 
extend the present procedure to the mean field equations 
governing a Boson- Fermion mixture of atomic gases [37| . 
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